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Abstract. A modified Chaplygin gas model (MCG), pmcg/pmcgo = [B s + (1 - B s )a" 3(1+B)(1+a) ] 1/(1+o) , 
as a unified dark matter model and dark energy model is constrained by using current available cosmic 
observational data points which include type la supernovae, baryon acoustic oscillation and the seventh 
year full WMAP data points. As a contrast to the consideration in the literatures, we do not separate 
the MCG into two components, i.e. dark mater and dark energy component, but we take it as a whole 
energy component-a unified dark sector. By using Markov Chain Monte Carlo method, a tight constraint 



is obtained: a = 0.000727 
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1 Introduction 

Under the assumption of Einstein gravity theory, from the 
observed energy components T°^f and the geometric struc- 
ture G tIlVl the un-observed energy component is defined as 
a dark fluid 



rpdark 
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rpobs 



(1) 



In the literatures, a dark fluid is usually separated into 
dark matter and dark energy parts, for recent reviews on 
dark energy please see [TJUJOIjmGiHHlG] • F° r gravitational 
attraction, the dark matter component is responsible for 
galaxy structure formation, and dark energy provides re- 
pulsive force for current accelerated expansion JSUH] • How- 
ever, we can not observe the dark sectors directly in cos- 
mology. From the above equation ([T]), one can easily see 
the degeneracies between dark fluid and gravity theory 
and can understand the reasons of why a kind of modified 
gravity theory, for example f(R) gravity theory, can re- 
alize current accelerated expansion. But what is the con- 
tent of the dark fluid? Be honest, one can give a long 
list of candidates which include dark matter and dark en- 
ergy models, but the last word is still empty until now. If 
the Higgs particles were founded in LHC, the scalar field, 
quintessence, would be a competitive candidate. But an- 
other possibility is that the dark fluid is a mixture of dark 
matter and dark energy components, or it is just one ex- 
otic unknown fluid. This property, dubbed dark degener- 
acy, has been discussed by the authors [10] . A virtue of the 



unified dark fluid model is that the so-called coincidence 
problem is removed completely. 

A modified Chaplygin gas (MCG) model [TT], which is 
a unified dark matter and dark energy model, is an exam- 
ple of dark degeneracy. MCG model is a variant of Gener- 
alized Chaplygin gas (GCG) [T3] which is a generalization 
of Chaplygin gas (CG) [13]. MCG has been discussed in 
many perspectives extensively [T^lfTSHTfj] . 

In our previous work [TS], the energy density of MCG is 
decomposed into two components: one part, which evolves 
with the calling relation a -3 , is dark matter; and the re- 
maining part is dark energy. With this decomposition, by 
using 182 Gold SN la, 3-year WMAP and SDSS BAO, 
the best fit values of the model parameters were obtained: 
B s = 0.822, a = 1.724 and B = -0.085. In our work [TO], 
we used CMB shift parameters, BAO, SN la Union 2, ob- 
servational Hubble data and cluster X-ray gas mass frac- 
tion (CBF) to constrain the model space: B s = 0.7788t 0736 
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B = 0.00189±g;™^ 3 (la) ±g:SooxS ( 2<T )> where MCG was 
taken as a unified dark sector. But in the CBF, BAO 
and CMB constraints, an effective matter density Q m = 
n b + (1 - O b - Q r - J? fc )(l - B,) 1 /* 1 ^' was used. So in 
some senses, the potential decomposition was employed 
too. When one just considers the background evolution 
and doesn't need any kind of definition of Q m , one doesn't 
worry about this issue. But one has to take into account 
decomposition carefully, when the perturbations of energy 
component is involved. The problem comes from the fact 
that the decomposition is not unique. One can give many 
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kinds of decomposition, because of the lack of physical 
principle to do a decomposition. The worst is that the 
evolutions of perturbations strongly depend on the de- 
composition and definition of dark matter. As a contrast, 
in this paper, we shall not take any decomposition like 
that in our previous work [T51ITC] . i.e. we take MCG as a 
whole energy component entirely. Furthermore, the per- 
turbations of MCG will also be included. In our previous 
works [15,16], the background evolution information and 
the shift parameters R and l a and but not the full 
information from CMB was used. One can expect that, 
by combining the full information from CMB, a tighter 
constraint would be obtained. So, in this paper, we shall 
constrain the background evolution by using BAO and SN 
la data points, and its perturbation evolution by using the 
full CMB data. Finally, the model parameter space is ob- 
tained. 

This paper is structured as follows. In section [21 the 
equation of state (EoS) and adiabatic sound speed of MCG 
are shown. Meanwhile, the Friedmann equation and per- 
turbation equation are given. In section [31 the constraint 
method and results are presented. A summary is given in 
Section g] 



2 Main equations in modified Chaplygin gas 
model 

The MCG is characterized by its equation of state (EoS) 



Pmcg = Bpmcg - A/ pi 



MCG 



(2) 



where B, A and a arc model parameters. It is obvious 
that GCG is recovered when the value of B is zero. And, 
a cosmological constant A is reduced when a = — 1 and 
A = 1 + B. Also, if A = is respected, MCG looks like a 
perfect fluid with EoS w = B, for example a quintessence 
model. However, if MCG as a unified dark sector, it would 
not happen. Because it would look like a combination of 
cold dark matter and a simple cosmological constant. In 
general, for a spatially non-flat FRW universe, the metric 
is written as 



ds 2 = -dt J 



a 2 (t) 



1 



1 — kr 2 



dr 2 
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( 3 ) 

where k = 0, ±1 is the three-dimensional curvature and a 
is the scale factor. Considering the energy conservation of 
MCG, one can rewrite the MCG energy density as 



equation 

H 2 = H 2 {n b a- 3 + Q r ar A + /2 fc a~ 2 

+ (1 - Q b - Q r - Q k ) \B S + (1 - B s )a-^ 1+B ^ 1+ ^ 



where H is the Hubble parameter with its current value 
H = 100/ikm s _1 Mpc _1 , and fli (i — b : r : k) are dimen- 
sionless energy parameters of baryon, radiation and effec- 
tive curvature density respectively. In this paper, we only 
consider the spatially flat FRW universe. 

To study the effects on CMB anisotropic power spec- 
trum, the perturbation evolution equations for MCG would 
be studied. We treat MCG as a unified dark fluid which 
interacts with the rest of matter purely through gravity. 
With assumption of pure adiabatic contribution to the 
perturbations, the speed of sound for MCG is 

c 2 = ^- = j j = -aw+{l + a)B, (6) 



where w is the EoS of MCG in the form of 

B K 



w = B - (1 + B) 



B s + (l-B s )a- 3 ( 1 + B )( 1 +«)' 



(7) 



In the synchronous gauge, using the conservation of energy- 
momentum tensor TJf = 0, one has the perturbation 
equations of density contrast and velocity divergence for 
MCG 

Smcg = -(1 + w)(0mcg + ^) - 3H{c 2 s - w)5 M cg (8) 

c 2 

Omcg = -H(l - 3c 2 s )9 M cg + -^k 2 6 MC G - k 2 a M M 

1 + w 

following the notation of Ma and Bertschinger [17]. For 
the perturbation theory in gauge ready formalism, please 
see |18j . The shear perturbation o~mcg — is assumed 
and the adiabatic initial conditions are adopted in our 
calculation. 

To keep the perturbation evolution stable, the positiv- 
ity of sound speed c 2 is demanded. Actually, to protect 
the causality, an upper bound c 2 < 1 is also needed. If 
one has the assumption of a > and < B s < 1, one can 
find easily that c 2 is non-negative if B respects to the in- 
equality B > —aB s /(l + aB s ). However, this assumption 
is unnatural. Here we shall take a as a free model param- 
eter and remove any unnatural assumption. The detailed 
methodology for keeping the positivity of sound speed c 2 
will be discribed in the next section. 



Pmcg — Pmcgo 



B s + (1- B s )a 



-3(l+B)(l+a) 



(4) 



3 Constraint method and results 



for B ^ -1, where B s = A/(l + B)p\jg G0 . We shall take 3.1 Method and data points 



B s , B and a as MCG model parameters in this paper. 
Form Eq. (0}, one can find that < B s < 1 is demanded 
to keep the positivity of energy density. If a — and B = 
in Eq. Q, the standard ylCDM model is recovered. Tak- 
ing MCG as a unified component, one has the Friedmann 



We perform the observational constraints on parameter 
space by using Markov Chain Monte Carlo (MCMC) method 
which is contained in a publicly available cosmoMC pack- 
age [12], including the CAMB [3D] code to calculate the 
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theoretical CMB power spectrum. We modified the code 
for the MCG as a unified fluid model with its pertur- 
bations included. The following 8-dimensional parameter 
space is adopted 

P = {io b , S , t, a, B, B s ,n s ,log[10 10 A s }} (10) 

where uib = fibh 2 is the physical baryon density, 0s (mul- 
tiplied by 100) is the ration of the sound horizon and angu- 
lar diameter distance, r is the optical depth, a, B and B s 
are three newly added model parameters related to MCG, 
n s is scalar spectral index, A s is the amplitude of of the 
initial power spectrum. Please notice that the current di- 
mensionless energy density of MCG flucc is a derived 
parameter in a spatially flat (k = 0) FRW universe. So, it 
is not included in the model parameter space P. The pivot 
scale of the initial scalar power spectrum k s o = 0.05Mpc _1 
is used. We take the following priors to model param- 
eters: uj b G [0.005,0.1], S G [0.5,10], r <E [0.01,0.8], 
a G [-0.1,0.1], B G [-0.1,0.1], B s G [0,1], n s G [0.5,1.5] 
and log[10 10 Aj] G [2.7,4]. In addition, the hard coded 
prior on the comic age lOGyr < to < 20Gyr is imposed. 
Also, the weak Gaussian prior on the physical baryon den- 
sity uib = 0.022 ± 0.002 [5T] from big bang nucleosynthesis 
and new Hubble constant Hq = 74.2 ± 3.6kms~ 1 Mpc~ 1 
[22 are adopted. 

As is seen in Eq. the expression of c 2 which con- 
tains model parameters a, B, B s and scale factor a, is 
complicated. Giving an explicit range of model parame- 
ters to keep c 2 nonnegative is really difficult. Maybe, in 
some senses, it is impossible. To circumvent the problem, 
we take the code as a black box, and hard code the condi- 
tion c 2 (a) > 0. It means that, in very sampling, the code 
checks whether the value of c 2 is negative or not at first. 
If it is negative, the parameters combination is thrown 
away and a new one is generated randomly. If it is not 
negative, the evolutions of background and perturbation 
for MCG will be calculated. In this way, the output pa- 
rameters space holds the condition c 2 > 0. And we do not 
worry about its upper bound, because the value of c 2 is 
very small and near to zero. So the check condition be- 
haves like a filter. We can check the evolution of c 2 with 
respect to scale factor a once the final result is obtained. 
As is shown in the following section, please see Figure [3] 
it really works. 

The total likelihood C oc e~ x / 2 is calculated to get 
the distribution, here % 2 is given as 



it is not true in our case. So, we find Zd numerically from 
the following integration [26] 



2 _ 2 
X — XCMB 



Xbao 



XSN- 



(11) 



The CMB data include temperature and polarization power 
spectrum from WMAP 7-year data [23 as dynamic con- 
straint. The geometric constraint comes from standard 
ruler BAO and standard candle SN la. For BAO, the 
values {r s (zd)/Dv (0.2), r s (z t ;)/Dy(0.5)} and their inverse 
covariant matrix [21] are used. To use the BAO informa- 
tion, one needs to know the sound horizon at the redshift 
of drag epoch Zd- Usually, Zd is obtained by using the 
accurate fitting formula [22] which is valid if the matter 
scalings pb oc a~ 3 and p c oc a -3 are respected. Obviously, 



TiVd) = / drj'id 

' Zd dr\x e (z)a T , 

dz- = 1 

u eta it 



(12) 



where R = 3pb/4:p 7 , ut is the Thomson cross-section and 
x e (z) is the fraction of free electrons. Then the sound hori- 
zon is 

r s {zd)= \ dr]c s (l + z). (13) 



where c s = + R) is the sound speed. We use the 

substitution [55] 



, r s (z d ) . . 
d z : -, J s(Zd), 
r s {Zd) 



(14) 



to obtain unbiased parameter and error estimates, where 
d z = r s (zd)/ Dy(z), f s is evaluated for the fiducial cos- 
mology of Ref. [24], and id is obtained by using the fit- 
ting formula [3S] for the fiducial cosmology. Here Dy(z) — 
[(1 + z) 2 D\cz/H(z)] 1 ^ 3 is the 'volume distance' with the 
angular diameter distance Da- The 557 Union2 data with 
systematic errors are also included [37] ■ F° r the- detailed 
description of SN, please see Refs. [28] . 



3.2 Fitting Results and Discussion 

We generate 8 independent chains in parallel and stop 
sampling by checking the worst e- values [the variance(mean) /mean(\ 
of 1/2 chains] R — 1 of the order 0.01. The calculated re- 
sults of the model parameters and derived parameters are 
shown in Table. [T] where the mean values with la and 
2<T regions from the combination WMAP+BAO+SN are 
listed. The minimum x 2 is 8004.630 which is smaller than 
that xL™ = 8009.116 for /1CDM model with the same 
data sets combination. Correspondingly, the contour plots 
are shown in Figure [1] 



From the Table Q] and Figured] one can clearly see that 
a tight constraint is obtained when the full information of 
CMB data is included. For the small values of a and B, 
one finds that MCG model is very close to ylCDM model. 
And the current data slightly favor MCG model. Using the 
mean values of model parameters, we plot the evolution 
of EoS w(a) and the speed of sound c 2 (a) for MCG with 
respect to scale factor a in Figure [5] From the left panel 
of Figure [U the MCG behaves like cold dark matte at 
early epoch. The right panel of Figure [H shows that the 
value of the sound speed c 2 (a) of MCG is small positive 
number and varies with scale factor. The small values of 
c 2 (a) make it possible to form large scale structures in our 
universe. Also, the positivity of c 2 (a) is really guaranteed 
by the 'filter' in sampling. 



4 Lixin Xu et al.: Modified Chaplygin Gas as a Unified Dark Matter and Dark Energy Model and Cosmic Constraints 



Prameters 



Mean with errors 



e 

T 

a 
B 

B s 
n a 

log[10 10 A 3 ] 
Omcg 
Age/Gyr 

n b 

Zre 

H 



— n n 99«+ l, 4)oob40+o.ooi08 — 

U.UZZO_ .0005360.00104 

i n/iQn+ - 00251 + - 00502 

L.UiyU_ 002521-0. 00494 

n ns«s+°- 00659 + - 0255 

U.UOOO_ .00735-0 .0238 

n nnn7O7+ - 00142 + - 00391 

U.UUU I Z I _o. 00140-0. 00234 

n nnn777+ ' 000201 + ' (l(,(l ' rl5 

U.UUU l I I -0.000302-0.000697 
fl 7QO+0.0163+0.0307 
u - loz -0. 0162-0. 0329 
fl Q.07+0. 01451+0. 0287 
U.HOI _ .0146-0 .0285 
o my . +0.03:«+0.0li95 
d-UOM_ 0i 0339-0. 0651 
fl 01=7+0. 00182 + 0. 00359 
U.yo / -0.00184-0.00363 
1 3 ROQ+0. 120+0. 228 
io - ozy -0. 118-0. 235 
11 fl |-; + 0. 00184+0. 00808 
U.UIO.o. 00182-0. 00359 
10 c 9zl +l. 2152+2. 417 
1U.OZ4_j .216-2.400 
79 +1 679+3.361 
' z ' oul -l. 690-3. 223 



Table 1. The mean values of model parameters with lcr and 
2cr errors from the combination WMAP+BAO+SN. 



0.0 

-0.2 

-0.4 

-0.6 

-0.8 
0.0 

0.0020 r 
0.0015 



0.2 



0.4 



0.6 



0.8 



1.0 



0.0010 



<m 


l/v 

i 


l/\ 






| 




IAI 


i\m 


|©> I 


!#> 


\^\ 


mi 


| 

l«> 1 


."[%> 

1 


:\ ^)| 
1^1 



Fig. 1. The ID marginalized distribution on individual param- 
eter and 2D contours with 68% C.L. and 95% C.L. by using 
CMB+BAO+SN data points. 



To understand the effects of model parameters to the 
CMB anisotropic power spectra, we plot the Figure |3l 
where one of three model parameters a, B and B s varies 
in the first three panels, and the other relevant parame- 
ters are fixed to their mean values as listed in Table [TJ The 
upper left and right panels of Figure [3] show the effect of 
parameter a and B to CMB power spectra respectively. 
The model parameters a and B modify the power law of 
the energy density of MCG, then they make the gravity 
potential evolution at late epoch of the universe. As re- 
sults, one can see Integrated SachsWolfe (ISW) effect on 
the large scale as shown in the upper left and right panels 
of Figure [3] In the early epoch, MCG behaves like cold 
dark matter with almost zero EoS and speed of sound c 2 s 
as shown in Figure [U therefore the variation of the val- 
ues of a and B will change the ratio of energy densities 
of the effective cold dark matter and baryons. One can 
read the corresponding effects from the variation of the 
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Fig. 2. The evolutions of EoS w(a) and sound speed (? s {a) for 
MCG with respect to scale factor a. 



first and the second peaks of CMB power spectra. The 
parameter B s is related with the dimensionless density 
parameter of effective cold dark matter i? c n,. Decreasing 
the values of B s , which is equivalent to increase the value 
of effective dimensionless energy density of cold dark mat- 
ter, will make the equality of matter and radiation earlier, 
therefore the sound horizon is decreased. As a result, the 
first peak is depressed. The lower right panel shows CMB 
power spectra with mean values listed in Table[T]for MCG 
and ylCDM model, where the black dots with error bars 
denote the observed data with their corresponding uncer- 
tainties from WMAP 7-year results, the red solid line is 
for MCG with mean values as shown in Table [TJ the blue 
dashed line is for ylCDM model with mean values for the 
same data points combination. And the green doted line is 
for ylCDM model with mean values taken from [29 with 
WMAP+BAO+i?o constraint results. One can see that 
MCG can match observational data points and /1CDM 
model well. This is the evidence of dark degeneracy. 



4 Summary 

In summary, we perform a global fitting on MCG model, 
which is treated as a unified dark matter and dark energy 
model, by using MCMC method with the combination of 
the full CMB, BAO and SN la data points. As a contrast 
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to the reports in the literatures, we take MCG as an entire 
energy component and without any decomposition. Tight 
constraint is obtained as shown in Table [T] and Figure 
[TJ The MCG model can match observational data points 
and ylCDM model well. This is the dark degeneracy. For 
the small values of model parameter a and B, one can 
conclude that MCG model is very close to ylCDM model. 
And the current data favor MCG model slightly. 
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Fig. 3. The CMB C ; power spectrum v.s. multiple moment 
I. The upper left, right and lower left panels show the effects 
of three model parameters a, B and B s to CMB tempera- 
ture anisotropic power spectra respectively, in each case the 
other relevant model parameters are fixed to the mean val- 
ues as listed in Table [T] The lower right panel shows CMB 
power spectra with mean values listed in Table[T]for MCG and 
ylCDM model, where the black dots with error bars denote 
the observed data with their corresponding uncertainties from 
WMAP 7-year results, the red solid line is for MCG with mean 
values as shown in Table [T] the blue dashed line is for /1CDM 
model with mean values for the same data points combination. 
And the green doted line is for /1CDM model with mean values 
*„1»™ f,.™, [onl „,;+u ixn\,f Ar> D»A i u „„„„+„„;„+ ,-„„,, i+„ 



